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Abstract 

We examine the quantisation of a collective Hamiltonian for the two-baryon system 
derived by us in a previous paper. We show that by increasing the sophistication 
of the approximations we can obtain a bound state - or a resonance - not too 
far removed from the threshold with the quantum numbers of the deuteron. The 
energy of this state is shown to depend very sensitively on the parameters of the 
model. Subsequently we construct part of a collective Hamiltonian for the three 
baryon system. Large-amplitude quantum fluctuations play an important role in the 
intrinsic wave function of the ground-state, changing its symmetry from octahedral 
to cubic. Apart from the tetrahedron describing the minimum of the potential, we 
identify a "doughnut" and a "pretzel" as the most important saddle points in the 
potential energy surface. We show that it is likely that inclusion of fluctuations 
through these saddle points lead to an energy close to the triton's value. 



1 Introduction 



In a recent paper [1] we have investigated a numerical procedure to determine 
a collective coordinate manifold for the Skyrme model, and have exhibited 
many of the parameters in the collective Hamiltonian. At almost the same time 
Leese, Manton and Schroers [2] published a paper studying the quantisation 
of the Skyrme model in the attractive channel. Even though their collective 
Hamiltonian was not determined self-consistently, they obtained some very 
interesting results. They found a bound state at 6 MeV binding energy, with 
the quantum numbers of the deuteron. More importantly they found a pion 
density remarkably similar to that in the deuteron. 
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In this paper we concentrate on a description of the deuteron using our collec- 
tive Hamiltonian. We have a lot of additional information, that can be used to 
construct several approximations of increasing sophistication to the dynamics 
of the Skyrme model in the 12-dimensional manifold thought to describe the 
dynamics of the deuteron. We shall also exhibit the dependence on the choice 
of parameters in the Skyrme model, and shall argue that the result of Manton 
and collaborators was accidental, and that the best one can hope for is an 
energy not too far from zero, maybe with an error of 30 MeV or so. 

After having analysed the deuteron in great detail we turn our attention to 
the triton and 3 He. A while ago, Carlson [3] has studied the application of 
the Skyrme model to the system of three baryons (see also Ref. [4] where the 
B = 3 system is studied using the approximate Yang-Mills instanton-induced 
(YMI-induced) form for the Skyrme fields). In these studies one starts from the 
minimum energy solution where the baryon density has tetrahedral symmetry 
(see Fig. la). One then makes the approximation that this minimum describes 
the triton. This approximation leads to a tremendous over-binding, which was 
attributed to the neglect of simple quantum effects, especially vibrational zero- 
point motion. It has also been mentioned [5] that anharmonic modes might 
play a role. 



Fig. 1. A plot of a surface of constant baryon density in (a) the tetrahedron , (b) 
the B = 3 doughnut and (c) the pretzel. 

In this paper we show that the low-energy potential "landscape" has a lot of 
structure, which has both important qualitative and quantitative effects. The 
most salient features in this landscape, apart from the tetrahedral solutions, 
are two meta-stable states at slightly higher energy, the B = 3 doughnut and 
the "pretzel" which has a planar symmetry similar to the doughnut but has 
two holes (see Figs, lb and c). 

The approach taken is this work is based on the techniques discussed in great 
detail in Ref. [1], where we show how to study large amplitude collective 
motion in the Skyrme model, using the YMI-induced forms of the fields for 
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simplicity. This is based on a mode-following approach, where we start from 
the harmonic fluctuations around a stable solution, and follow those into the 
anharmonic regime. One advantage of our approach is that such anharmonic 
modes always run through the extrema of the potential energy. The main 
limitation of the current work is the use of the instanton induced form of the 
Skyrme fields, see also Ref. [4]. 

The paper is organised as follows. In Sec. 2 we introduce the Skyrme model, 
and discuss the instanton-induced approximation to the dynamics of the model. 
We also discuss the relation of the scaled Skyrme units to standard units. We 
then discuss, in Sec. 3, how the study of harmonic fluctuations can be used 
to calculate improved observables. In the next section, Sec. 4, we succinctly 
recapitulate the properties of the B = 1 hedgehog within the YMI-induced 
Ansatz. We then discuss the quantisation of the B = 2 system in Sec. 5. In 
Sec. 6 we discuss the structure of the potential landscape, and how the different 
structures are connected. In the next section (Sec. 7) we try to estimate what 
the effects of the structure of the potential landscape are on the ground-state 
of the triton. Finally, we draw some conclusions in Sec. 8. 



2 Skyrme Lagrangian, Skyrme units, and choice of parameters 



2.1 The model 



The "standard" Skyrme model is based on the non-linear sigma model, ex- 
tended by a quartic interaction term and a pion mass term. The instanton- 
induced form [6,7] of the Skyrme fields can only be used in the (chiral) limit 
of zero pion mass, where the model is defined by the Lagrange density 

f 2 1 

£ = J -±Tr{d ll Ud> l rf) + -Tr[U^d^U, d v U\[U ] d^U , U^d u U], (1) 

4r O^O 



where U is a unitary two-by-two matrix-valued field satisfying the bound- 
ary condition U — 1 at infinity. As has been discussed many times before, 
this model has a topologically conserved quantum current. The charge of this 
current is identified with baryon number B. 

If we rescale the units of time and length, x — > xj (ef n ), (the so-called Skyrme 
units), the Lagrange density takes on the slightly more convenient form 

£ = k f-Tr(d^Ud fl U r ) + — Tr[U^d^U, U ] 'd v U][U ] d»V ", U ] d v U]\ , (2) 
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where f n /e is the Skyrme unit of energy. 

Finally the Skyrme Lagrangian can easily be reformulated in terms of the 
Sugawara variables (Lie-algebra valued currents) L, 

L fl = U~ l d^U = iiy a , (3) 



and we have 



£ = 



1 jajn , 1 

2 V" ^ 4 



i°W*i vb 



(4) 



2.2 Instanton-induced Skyrme fields 



As discussed in [7] one can derive a Skyrme field from a Yang Mills instanton 
field (called YMI-induced in this paper) by integrating the time component of 
the gauge potential, 



U(x) = CS { Pexp 



-A 4 (x, t)dt 



CK 



(5) 



Here S is a constant matrix, chosen such that U decays to 1 at infinity, and 
C describes an overall grooming. For the current work, where we shall only 
consider C near the identity, it is convenient to parametrise 



— * 

C = exp(ir • 9). 



(6) 



For the Jackiw-Nohl-Rebbi (JNR) instanton of charge k we have [8] 



A 4 (x,t) = --f-T, P=E|„. v 
Z p l=1 \X — Aj| 



(7) 



and we should use S = —I, to obtain a field of positive baryon number k. 

In order to solve for the YMI-induced value of U we convert the integral (5) 
to the solution of a differential equation. First introduce 

U(x,t) = CS jpexp j -A 4 (x,t)dt^ CK (8) 
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This function satisfies the differential equation 



d T U(x,r) = —U(x,t)A 4 (x,t), 



(9) 



with initial condition U(x,—oo) = S. The function U(x) is obtained as the 
limit for r — > oo of U. We can work directly with the Sugawara variables L^, 
Eq. (3). These can be calculated as the large r limit of a quantity L M , which 
satisfies the differential equation 



drL^x, t) = [A 4 (x, t), Lh(x, t)] - d^A 4 (x, r). 



(10) 



Here the boundary condition is L^(x, — oo) = 0. By differentiating the dif- 
ferential equations (10) we can obtain expressions for derivatives of L M with 
respect to the instanton parameters A/ and X h which will be needed later. 

The field U defined in Eq. (5) does not have an explicit time dependence. 
There is an implicit dependence, due to a possible variation of the instanton 
parameters as well as the unitary matrix C (parametrised by the three angles 

— * — * 

9) with time. Let us denote the parameters {\ h X h 9} collectively by £. We 
then have 



L = i a U ] d^U = i a L : , 



(11) 



If we substitute this in the Lagrangian we obtain the form 

L = T - V, 



(12) 



with 



(13) 



and 



B <*p = E K<*Kp + 2 



E'»E'#-E(ETOE# 6 



Lb 



(14) 
(15) 



The Lagrangian is quadratic in the time-derivatives due to the special nature of 
the Skyrme model. This also means that we have broken the Lorentz invariance 
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of the original equations, so that (11) can only be used to describe adiabatic 
(small velocity) motion. 

In the YMI-induced Ansatz we have thus replaced the general matrix U (x) , 
with an infinite number of parameters, by a form parametrised by 3 + 5(k + 1) 
parameters, CU(x\£)C^. 

2. 3 Skyrme units and model parameters 
Table 1 

Different parameter sets as used in the Skyrme model 

set f w e 

ANW 64.5 MeV 5.45 

ANW 81.7 MeV 6.95 

R 90 MeV 4 

We take three reasonable parameter sets in order to be able to compare the 
effect of changes in the parameters. The standard values for zero pion mass 
are those fitted by Adkins, Nappi and Witten [9] to the N and A masses, and 
are given in table 1 under the label "ANW". Since in the instanton-induced 
Skyrmion the moment of inertia differs considerably from its exact value, we 
propose a different set of parameters under the label "ANW", that give the 
exact N and A masses for the instanton-induced fields. Finally we use a third 
set of parameters with realistic /„- and e. Of course we now no longer reproduce 
the masses. 
Table 2 

Our Skyrme units and their values 

parameter/ units SU ANW ANW R 

[length] l/(/w e ) 0.562 fm 0.347 fm 0.548 fm 

[energy] f v /e 11.8 MeV 9.28 MeV 22.5 MeV 

h e 2 29.7 48.3 16 

Now it may be useful to recapitulate our units, which we have done in table 
2. Note that H is not one in Skyrme units, but c is. 



3 Harmonic expansions 

In these notes we shall use a harmonic expansion around (quasi-) stable solu- 
tions of the Skyrme model. Starting from the Lagrangian (1) we can use the 
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instanton-induced fields to formulate a harmonic approximation of the form 
(we use £ for the deviation of the coordinates from their equilibrium value) 

Hho = E + ^B a p£ a £p + \V a p£ a £". (16) 

We can diagonalise the harmonic Hamiltonian by making a linear transforma- 
tion to new coordinates (see e.g. Ref. [10]) 

(17) 

where / and g are the left and right eigenvectors of the matrix 

M a p = B^V n p, (18) 

obeying the normalisation condition 

gX = (19) 

and with eigenvalues (hu^) 2 . The Hamiltonian can then be cast in the diagonal 
form 

H = E + ]T 'tiw^n, + ]T K'pI (20) 

There the first sum runs over all non-zero eigenvalues, and the second one only 
over the zero modes. Of course the quantised form of this same Hamiltonian 
can easily be obtained by replacing n M by n M + 1/2, interpreting n M and as 
operators. As is well known, see e.g. Ref. [11], there is a correction to E (due 
to normal ordering of the excitation operators and the non-zero expectation 
value of the operators pi in the ground state), that usually more than cancels 
the harmonic zero-point energy. Unfortunately, these contributions are not 
easily evaluated in a field theoretical context. See, however, Ref. [13]. 

Using similar techniques, one can easily obtain expressions for the ground state 
expectation values of observables. Of more direct interest when interpreting 
the harmonic modes is the change in observables related to a given mode. We 
shall only discuss quantities independent of the velocities. We expand 

o = o +eo, a + ie^o, a] p, (2i) 
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O jCt etc. are nothing more than (covariant) derivatives w.r.t. £ Q . Of course we 
can apply the chain rule to convert these to derivatives w.r.t. g M . To linear 
order one finds, for small 5£, 

50 = 5CO, a = 6<fglO, a . (22) 



From which we see that 

O,, = g a »0, a (23) 



is a good definition of the change of O associated with the mode //. This 
expression is very useful in interpretation of the harmonic modes. We can use 
it to define the change in baryon number, which is highly enlightening. 

In order to calculate the expectation value of the operator O in the ground 
state we now quantise the coordinate £ (or rather q) using the standard har- 
monic oscillator rules 



Q 



^ = -\¥) ~V^' (24) 

where we have used the fact that B a p and V a p can be diagonalised simultane- 
ously (with eigenvalues b and v, respectively). We have hu^ = (b^v^) 1 / 2 . 

The quantised form of the observable O now becomes 
1 /b^\ 1/2 

O = O + -J2'[-j fZfHO, a .f, + normal ordered terms. (25) 



The second term can now interpreted as the zero-point- mot ion contribution to 
O. Notice that everything above was done in the intrinsic frame (also called 
the body-fixed frame). In order to obtain the result for an observable in the lab 
frame it still has to be averaged over with a suitable rotational wave function 
[11]. 
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4 Properties of the B = 1 hedgehog 



The B = 1 collective Hamiltonian is given by zero-modes only, 

w = B ° + !<' 2+s2 > + 2s p2 - < 26 » 



The parameters in if are given in table 3, where we both show the results 
for the YMI-induced approximation, and for the numerical solution of the full 
problem. 

Table 3 

Properties of the B = 1 system. 





instanton-induced exact 




73.6 f v /e 


73.0 fje 


h 2 /A 


1/140.1 Ue 3 


1/106.6 f n e 3 



5 Quantisation of the B = 2 system. 



In this section we shall use the results obtained in Ref. [1] for the B = 2 system 
to make a set of approximations to the energy of the deuteron of various levels 
of sophistication. We shall start with a naive harmonic analysis, and end with 
an approximation that goes beyond the calculation by Leese et al [2]. 



5.1 Harmonic analysis. 



In this section we shall discuss the harmonic analysis around the doughnut 
in what has been called M.yi by Manton et al. This is a manifold that in 
the asymptotic regime has just enough collective coordinates (12 = 6x2) 
to be able to project the individual Skyrmions onto nucleons. Three of those 
describe the centre-of-mass motion, and do not couple. The remaining 9 modes 
describe the intrinsic dynamics. 

The issue is the understanding of bound states in this channel, where we 
restrict ourselves to the harmonic approximation. Since the doughnut has a 
continuous symmetry (generated by = 21^ + Jg nt ) we typically find 

two-dimensional pairs of modes when the quantum number m™* is changed, 
and ID when it remains zero. For modes breaking the symmetry we find a 
centrifugal term of the form h 2 / ex 2 (Mg nt ) 2 (where x is a coordinate describing 
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the deviation from equilibrium). This leads to a 2D harmonic oscillator type 
spectrum, which can be given in the following form 



E = E + ]T hwifa + mf + 1)+]T hu t (n t + 1/2) 



2D ID 

+|-(/ 2 - (If) 2 ) + ^(J 2 - {J'Tf) + ^(I'f - 24 nt ) 2 - (27) 

Let me first discuss the classification of the modes: 
Table 4 

Harmonic frequencies around the B = 2 solution in the YMI-induced form. 
Classification deg. |Am3| Uuj /(/^e) 

translational zero modes (x — y) 2 1 

rotational zero modes (x — y) 2 1 

isorotational zero modes (x — y) 2 1 

translational zero mode (z) 10 

mixed iso/rotational zero mode (z) 1 

quadrupole mode 2 2 0.352 

dipole mode 2 1 0.419 

breathing mode 1 0.524 

highest modes 2 1 (?) 0.756 

Table 5 

Some properties of the B = 2 solution in the YMI-induced form. 

E 141.18 U/e 

h 2 /lj 1/114.2 f w e 3 

h 2 /lj 1/193.6 f n e 3 

h 2 /l c 1/65.2 f^e 3 

We can use these number when estimating the energy in M 12 relative to 2 
B = 1 Skyrmions. The ground state has intrinsic quantum numbers J = 
1, 1 = 0, J J™* = 0, If = 0. We find the following energy balance: 



E E, 



rot 



E B=2 = 141.2/,/e + 0.771/^ + ^f n e 



3 



2E B=1 = 147.2/Je + + f ^/.e 3 

AE = -6.0/ w /e + 0.771/^e + -.00819 /^e 3 



(28) 
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For our three sets of parameters this takes on the values 117.7, 117.0 and 88.2 
MeV. 

These estimates are all unbound, but let us analyse the situation a little bit. 
We know first of all that the moments of inertia are usually over-estimated 
by the YMI-induced Skyrme-fields. Suppose we reduce all moments of inertia 
by a factor 106/140, corresponding to the ratio of the exact over the YMI- 
induced value. The ground-state energy would then decrease by 27, 57 or 15 
MeV (the middle value should be disregarded, since this parameter set was 
already corrected for this effect). 

Furthermore we have ignored the fact that we know that the behaviour in 
the attractive channel direction (lowest 2 modes) is a lot less harmonic than 
in the other direction. The difference is due to the fact that in the attractive 
manifold the doughnut comes apart into two hedgehogs, and the doughnut is 
not strongly bound. The depth of the well is 6 S.U. which takes on the values 
71, 55.7 and 135 MeV for each of our three parameter sets; the harmonic 
energy quantum is 0.352e 2 S.U (the e 2 is really h in Skyrme units), which 
takes values 124, 200 and 127 MeV. In the other direction we expect to reach 
states with energies of the order of the repulsive channel (maybe up to 1 GeV), 
and the harmonic approximation should be good. 



148 




140 1 ■ 1 ■ 1 ■ 1 

10 20 30 

Q 

Fig. 2. The potential energy as calculated in the attractive channel (solid line). The 
dashed line shows the harmonic approximation to this potential, and the dotted line 
a Poschl-Teller fit 

In order to see what might be the possible effect of the limited well-depth 
we also calculate the binding energy of the lowest state by approximating the 
potential by a Poschl-Teller form [14] . This has the advantage that the problem 
is exactly soluble, and this potential has only two parameters, which can be 
chosen to be the depth and the harmonic frequency. If more information is 
known one of the generalised potentials discussed by Ginocchio [15] might be 
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of some use. For the time being we use the form 



V(x) = —v{y + 1)/ cosh 2 (o;x). 



(29) 



Here we have chosen our coordinates such that h 2 /m = 1. The parameters v 
and a can then be connected to the harmonic frequency and the depth of the 
potential 



z/ = -(-l + v /l + 8\/ /« 2 ) ) 
a = huo/ JV . 



(30) 



For this exactly solvable potential we then find that the ground state occurs 
at a zero-point energy 



En 



ftui 



Tiuj 



1 + 16 (Vo/hcu) 2 - 1 



(31) 



above the minimum of the potential. It is easy to see that this reduces to the 
harmonic oscillator value in the limit V ^> Tiuj. 

The reduction in the contribution for our lowest two modes would be 0.655, 
0.518 and 0.768 for our three parameter sets. This corresponds to 43, 96 and 30 
MeV reduction in binding. Taking all corrections together (we do not correct 
the middle value for the wrong value of the moment of inertia) we find 47, 
21 and 43 MeV. These values are reasonably close to zero, which is about the 
best one could hope for such a crude calculation. We have summarised our 
calculations in table 6. 
Table 6 

Summary of harmonic and related approximations for the deuteron energy relative 
to the B = 2 threshold, as discussed in the text. 





ANW 


ANW 


R 




harmonic 


118 MeV 


117 MeV 


88 


MeV 


corrected moment of inertia 


98 MeV 


117 MeV 


73 


MeV 


include anharmonicity 


47 MeV 


21 MeV 


43 


MeV 



5.2 Quantisation in M7 and extensions 



Leese, Manton and Schroers have constructed a Hamiltonian in the attractive 
channel, where one has 7 degrees of freedom. A quantisation of this Hamilto- 
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nian in the intrinsic frame led to a binding energy of 6 MeV. From our calcula- 
tion a similar Hamiltonian is available, and it would be nice to independently 
check these predictions, and we have all the necessary ingredient s.Q We nor- 
malise our collective coordinate such that the mass is 1 in Skyrme units, and 
we then have to quantise the intrinsic collective Hamiltonian (see also Ref. [2]) 



P 



(J. 



int\2 



int\2 



+ 



J: 



hit 



w(Q)l 



hit 



V 



a(Q) 2 KQ) 2 

ml ,2 ^jint^2 ( T \rA,\2 



(I u , 



A{Qf B{Qf C{Q)\ 



c(Q) 2 
+ V(Q). 



(32) 



Using the standard geometrical quantisation, and the intrinsic quantum num- 
bers compatible with the deuteron, I = 0, J = 1, J^ nt = 0, we obtain a "radial" 
equation of the form 



$(Q)=Ei>(Q). (33) 



Here g 1 ' 2 is the square-root of the determinant of the metric (mass matrix), 

g 1/2 = abcABC. (34) 



CD 
CD 

CD 

o 

H — ' 

CO 



0.003 



0.002 



0.001 




0.000 



10 



20 



30 



Q 



Fig. 3. The contribution of the rotational kinetic energy to the deuteron bound-state 
problem. 



2 Note that there is an error in the collective coordinate Q as presented in Ref. [1], 
as well as a factor of two difference in the inertial parameters. The results presented 
here use the corrected Q. 
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50 




Q 

Fig. 4. The metric factor g 1 / 2 in the kinetic energy, rescaled by a constant. 

The contribution of the rotational kinetic energy, which has units f n e 3 , is 
plotted in Fig. 3. We have subtracted the asymptotic behaviour - which cor- 
responds to the system separating into two hedgehogs locked in the attractive 
configuration. This is not the same as subtracting the kinetic energy of two 
hedgehogs in nucleon states. This would require that one includes the har- 
monic motion changing to zero- mode motion as the hedgehogs pull apart. 
This can (and should) only be done for scattering states [16]. For a bound- 
state calculations this part is missing (see below), one should only subtract so 
much rotational zero-point energy as to put the ionisation barrier at 0. 

The metric factor g 1 / 2 (scaled by a constant) is plotted in Fig. 4. It is very well 
fitted by a form aQ + f3Q 2 . The proportionality to Q for small Q agrees with 
the two-dimensional harmonic nature of the motion there. The asymptotic be- 
haviour agrees with the result by Leese et al. Mathematically one can interpret 
this as a system that changes it nature from two-dimensional at small Q to 3D 
at large separations. We have fitted all the relevant quantities as a function 
of separation, taking into account the exact behaviour for large separation. 
This allows one to solve the Schrodinger equation numerically. Some selected 
results are presented in table 7. 

In order to mimic the effects of quantising in the full 12-dimensional configura- 
tion space corresponding to full translational and iso-rotational zero modes of 
the two hedgehogs that appear at large separation we add the harmonic zero- 
point energies for these modes. Let us first study a correlation diagram, that 
shows all the non-zero modes as they develop from the doughnut, see Fig. 5. 
As one can see over there counting modes causes some problems. The lowest 
two modes correspond to the attractive channel, and the next two modes are 
the ones describing the two remaining fluctuations in M^i- The problem lies 
in the level crossing at Q — 8.8. If this were an avoided crossing, we would 
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Table 7 

The energy of the lowest bound state in units of f n /e as a function of e. The three 
columns denote the pure potential (no rotational energy), potential plus rotational 
energy and the additional effect of also adding the harmonic zero-point energies for 
the other modes in M\i- A dash indicates no solution. 

e potential rotational all 



1 


-5.67 


-5.67 


-5.49 


2 


-4.94 


-4.63 


-4.43 


3 


-3.95 


-3.83 


-2.47 


4 


-2.67 


-2.32 


-0.30 


5 


-1.28 


-0.65 




6 


-.0044 







0.8 



2 

CO 




-0.2 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

5 10 15 20 25 

Q 

Fig. 5. A correlation diagram, showing how the local harmonic frequencies are cor- 
related as we move from the doughnut in the attractive channel. 

happily follow the lowest mode. Unfortunately it is not: there is quite a lot 
of symmetry in the attractive channel and the two modes appear to have dif- 
ferent symmetries. This may be an artifact of the YMI-induced Ansatz, but 
it seems more likely that this is a real effect. Since this symmetry is closely 
linked to the attractive channel, it is highly probable that if we deform the 
system a little bit away from the attractive channel, we get an avoided cross- 
ing. The crossing point would then correspond to what has been termed a 
diabolic point in the theory of Berry's phases, and lead to a monopole type 
geometric phase. For the time being we ignore this effect and shall sum the 
zero-point energy of the two lowest modes, as shown in Fig. 6. 

Results including this harmonic contribution are give in the third column 
of table 7. Taking all the results together, we are forced to conclude that 
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0.20 
0.15 

o 

O 0.10 

E 

JC 0.05 
0.00 

10 20 30 

Q 

Fig. 6. The harmonic energy as added to the collective Hamiltonian. 

is almost impossible to find a bound state for a reasonable set of parameters 
when disregarding the harmonic motion, whereas when including the harmonic 
motion a value for e of at most four is allowed. This is a contradiction of the 
claim of 6 MeV binding energy by Leese et al [2] . They used a set of Skyrme 
parameters appropriate to a finite pion mass, however, without examining 
the parameter dependence. The fact that there is a bound state close to zero 
energy, or a resonance slightly above is the most important issue. One cannot 
consider the model to be right on scales of the order of a few MeV out of 
2GeV, especially not without including loops and the related counterterms 
[13]. 



6 The structure of the B = 3 collective manifold 

6.1 Harmonic modes around the tetrahedron 

For the B = 3 case we use the YMI-induced+JNR form as used by Leese 
and Manton [4]. This is not based on the most general instanton field but is 
probably sufficiently general for our purposes. 

The instanton field is generated by four poles, arranged in an tetrahedron. All 
weights are equal. We choose 

(A 1 ,X 1 ) = (l/4,-i?, -R, -R, 0), 
(\ 2 ,X 2 ) = (1/A,-R,R,R, 0), 
(A 3 ,X 3 ) = (1/4, J R,- J R, J R,0), 
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(35) 



where R is a parameter determining the size of the tetrahedron. 

Since all continuous symmetries are broken, we find 9 zero modes. The remain- 
ing modes can be classified as follows (The interpretation and classification of 
these modes is discussed in more detail in the appendix): 

10-11 hu = 0.318/ w e. 

These two modes correspond to the motion towards the B = 3 doughnut, 

See Fig. 7 for a sketch. 
12-14 hu = 0.332/^e. 

These modes are the beginning of a trajectory towards the pretzel mentioned 

in the introduction. 
15-17 huj = 0.469/ w e. 

It is highly probable that these modes describe the separation in a B = 2 

doughnut and a single hedgehog. 
18 hu = 0.483/^e. 

This is the breathing mode. 
19-21 hu = 0.813/„-e. These are the only modes that mix with iso-spin rotations. 

They describe a complex high-energy excitation. 

If we assume that we need 18 modes to describe the B = 3 collective coordi- 
nates, we notice that there is a clear gap in the spectrum just at that critical 
point. We have not used the most general instantons here, and it is unlikely 
but not impossible that this situation might change if we switch to the AHDM 
instanton [21]. 

Table 8 

Some properties of the B = 3 solution in the YMI-induced form. 



As said before, the lowest non-zero modes of the tetrahedron correspond to 
motion towards the B = 3 doughnut, where the line connecting any of the 
corners of the tetrahedron moves closer to the line connecting the other two 



Eq 

h 2 /i! 
h 2 /ij 
h 2 /i c 



206.6 U/e 
(113)- 1 Ue 3 
(357)- 1 f w e 3 
(-633)- 1 f n e 3 
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corners. When these lines cross we have reached the B = 3 doughnut which 
in our Ansatz is described by four poles arranged in a square. Actually one 
would expect three such modes, one for each coordinate axis. However, the 
sum of these modes is the breathing mode, which is in a different represen- 
tation of the tetrahedral group, and thus we find only two modes (see also 
our model discussion below). Following these modes naturally leads one to 
study the doughnut as well. It sits on the top of a saddle-point in energy. It 
has two unstable modes, whereas we expected only one corresponding to the 
path towards the tetrahedron. Indeed this is the most unstable mode. The 
other unstable mode leads to yet another saddle point. This solution has the 
double-hole structure of a pretzel. The potential energy of all these states, i.e., 
ignoring contributions from the rotational kinetic energy, is not very different 

-^tetrahedron = 206.56 f w /e, -E prctz ei = 211.54/^/e, doughnut = 214.06/ 7r /e, 
3E B=l = 220.8/Je, E B=2 + E B=1 = 214.8/ 7r /e. (36) 

As one can see from table 9 the sizes of the different configurations are differ- 
ent. Here we have given the values of the integrals 

(r?) = J r 2 B{r)d\ (37) 

which is the standard definition for these quantities. 
Table 9 

Sizes of the different shapes for B = 3 in S.U. 



solution 


(xx) 


(yy) 


(zz) 


(r 2 ) 


tetrahedron 


3.043 


3.043 


3.043 


9.129 


pretzel 


9.225 


2.436 


.916 


12.577 


doughnut 


4.974 


4.974 


.878 


10.836 



Let us first consider the paths connecting the different solutions. As we stated 
there is a path connecting the two independent tetrahedral configurations 
through a saddle point that looks like a doughnut. Actually there are three 
such paths, corresponding to the three "independent" doughnuts with rota- 
tional symmetry around each of the three coordinate axes. Such a path is very 
similar to the geodesic with twisted symmetry found in the scattering of three 
collinear monopoles by Houghton and Sutcliffe [17]. It corresponds to motion 
where the Skyrmion breaks apart into three B = 1 hedgehogs in an attrac- 
tive (and therefore twisted) configuration. The only thing is that in our case 
the Skyrmion does not go into into this form, but rather breaks up into two 
B = 3/2 "bananas" as shown in Fig. 8. Since solutions with fractional baryon 
number do not exist in the Skyrme model - they have infinite energy - this 
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Fig. 8. An example of the baryon density found when following the path through 
the doughnut and the tetrahedron, to a point beyond the tetrahedron 

leads to an infinite energy in the limit of infinite separation. It is plausible that 
the use of the JNR instanton, instead of the full AHDM data is responsible for 
this effect - It is hard to imagine how to describe the desired configuration in 
terms of the JNR data. For the motion in between the doughnut and hedgehog 
the JNR form appears to be reasonable. 

In Fig. 9 we show how the local harmonic frequencies change as we move 
from the doughnut to the tetrahedron and beyond. As usual we have fixed 
the normalisation of the coordinate Q by requiring that the mass (which we 
usually call B) equals one. There are no major exciting structures along the 
path, and one can see that the mode we followed (denoted by the triangles) 
interpolates smoothly between the tetrahedron and doughnut. Furthermore we 
find that several modes are doubly degenerate for all Q due to the symmetry 
in the configurations along the path. In this case ignoring the weights in the 
JNR instanton is no limitation; they are constant and all equal for all points 
along the path. The potential energy along the path is given in Fig. 10. 

Once one has the doughnut, one can study its unstable modes. It has two, 
the one discussed above, and one leading to the pretzel. This latter solution is 
described by a rhombic arrangement of the poles, so that for each doughnut 
we have two neighbouring pretzels, obtained by compressing any of the two 
diagonals of the square arrangements of poles describing the doughnut. In this 
case the weights in the JNR potential are also no longer equal, but this need 
not concern us here. The pretzel has one unstable mode, which connects to 
the tetrahedron, as we discovered through our mode following approach. 

A very important property of the path between the tetrahedron and the pretzel 
is that at all points the Skyrme fields have a trivial reflection symmetry under 
coordinate reflections, namely <f>(Rix) = Ri<p(x). 
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0.2 
0.0 
-0.2 

20 40 60 

Q 

Fig. 9. A correlation diagram, showing how the local harmonic frequencies are cor- 
related as we move from the doughnut to the tetrahedron (at Q « 20) and beyond. 
The lines with pluses and circles are doubly degenerate. The triangles are used 
to label the mode actually followed. The dashed line shows the new (zero-)mode 
appearing when we move away from the doughnut. Note that this connects to a 
non-zero mode at the doughnut. This is probably another deficiency of the JNR 
parametrisation. 

240 
230 

Z> 

CO 220 
> 

210 
200 

20 40 60 

Q 

Fig. 10. The potential energy as a function of the coordinate Q, as we move from 
the doughnut to the tetrahedron (at Q ~ 20) and beyond. Note the steep rise of the 
potential, above the "ionisation" limit of 220.8 S.U., when the system splits into 
two B = 3/2 "bananas". 

In Fig. 11 we show how the local harmonic frequencies change as we move 
from the pretzel to the tetrahedron. Actually we find that at the tetrahedron 
the mode followed (again denoted by squares) connects to the lowest triplet 
of modes in the doughnut. Since each of the modes of the triplet describes a 
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Fig. 11. A correlation diagram, showing how the local harmonic frequencies are 
correlated as we move from the pretzel (Q = 0) to the tetrahedron (at Q « 25). 

circuit through two pretzels and another tetrahedron, as denoted schematically 
in Fig. 12. Note that the lowest mode for the tetrahedron is the one connecting 



tetrahedron 



pretzel 



doughnut \ tetrahedron 



i pretzel 



Fig. 12. A schematic sketch of a few of the connections through a given doughnut. 
For each doughnut in one of the three coordinate planes, there are two associated 
pretzels, elongated along one of the two coordinate axes in the plane, and connected 
by a valley. There is a connection of the doughnut to two tetrahedra with different 
orientations. There also exists a closed valley connecting the tetrahedra through the 
pretzels. 



to the doughnut, which is the higher of the two saddle points. Obviously this 
figure only sketches one out of three connections, the whole structure of the 
collective manifold being highly complex. 



7 Quantisation of the collective coordinates 
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7.1 Harmonic approximation 



As in the case of the B = 2 system, the most straightforward way to get 
a "quantum" result for B = 3 is to quantise in harmonic approximation. 
Near the tetrahedron the collective Hamiltonian leads to an expression for the 
energy of the following form 



E = E + J2 M^ + l/2) 

modes 

h 2 h 2 h 2 -> - 

+ -—I 2 + ——J 2 + ——I ■ J. (38) 

2Xj 21 j 2Xc 



The ground state has intrinsic quantum numbers / = 1/2, J = 1/2, K = 0. 
Replacing the inner product of I and J with the square of the grand-spin 
operator K, we find the following energy balance: 

E ^^2huj E rot 

E B=3 = 206.6/^/e + 1.76/ ff e + .00353/^e 3 

(39) 

- 3E B+l = 220.SU/e + + | ^^e 3 _ 

AE = -14.2/ w /e + 1.76/^e - 0.0125/^e 3 



For our three sets of parameters this takes on the values 319.2, 386.5 and 
242.1 MeV. Here we once again see the effect of the over-estimate of quantum 
corrections. Notice that without these corrections the value would be —168, 
— 132 and —320 MeV, respectively. Since we expect anharmonicities to play a 
large role, we shall now study the effect of large amplitude fluctuations. 



7.2 Approximate large-amplitude dynamics 



As a fast and relatively simple way to study the modes one can assume that 
one can independently quantise the motion along the straight line and the 
ellipse, since these represent orthogonal modes at the tetrahedron. 

We first look at the paths through the doughnut. We have constructed a sim- 
ple model, based on the VSEPR model [12], as has been used in the chemistry 
of molecules. It is based on putting particles on a sphere, and letting them 
interact through a pairwise repulsive interaction, which is taken to be some 
power of the inverse of the separation of the particles. If we do that in our 
problem we should constrain the JNR poles on a sphere. The minimal energy 
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configuration will be a tetrahedron, and there will be a saddle point at the a 
square configuration of poles, corresponding to the doughnut. The configura- 
tion space is an eight-dimensional manifold. Three of those are the rotational 
zero modes. At the tetrahedron the remaining modes are a pair of degenerate 
modes describing the motion towards the squares and three further-non-zero 
modes. These last three are of no interest, but can be eliminated by impos- 
ing the additional constraint that the distance between two pairs of points 
is equal, and that the line connecting the two mid-points of each pair goes 
through the centre of the sphere. 




Fig. 13. Coordinates used in our VPSER-based model. 
With these constraints, using the coordinates specified in Fig. 13, we obtain 



fi = RD(Q) (sin 9 cos (ft, sin 9 sin 0, cos 9), (40) 

r 2 = RD(Q)( — sin 9 cos0, — sin 9 sin 0, cos 9), (41) 

r*3 = RD (Q) (sin 9 cos (p, — sin 9 sin 0, — cos 9), (42) 

r 4 = RD(Q)(— sin 9 cos0, sin 9 sin0, — cos 9). (43) 

(Here D(Q) is a rigid body- rotation corresponding to the three zero-modes). 
From this we can calculate the classical Lagrangian, 



i i<j 

= 2R 2 (Q 2 + 9 2 + sin 2 9$) 

-2/4" (l/(sin 2 9) n + l/(cos 2 9 + sin 2 9 sin 2 0) ra + 

l/(cos 2 # + sin 2 #cos 2 0)™) . (44) 
A lot of insight can be gained by plotting the potential as a function of 9 and 

— 

0. The kinetic energy is proportional to L , which allows the interpretation 
of the configuration space as a sphere. In order to plot the potential on this 
sphere, we have colour-coded the values of V. In Fig. 14 one can see this 
potential. We have also indicated the tunnelling paths from one minimum to 
the next. Here one can see why we have only two independent modes, and not 
three: one can not have three orthogonal paths in two dimensions. 
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Fig. 14. The potential energy (44) plotted as a function of and 4>. The white areas 
denote lowest potential energy, the black areas are where the potential diverges. 
The black lines show the tunnelling paths between different minima; the gray lines 
are the continuation of the black ones into the repulsive regime. 




Fig. 15. Orientation of the tetrahedra corresponding to each of the vertices of the 
cube in Fig. 16 

One has to be somewhat careful how to interpret this sphere. In essence 
there are only two independent orientations of the tetrahedron (the mini- 
mum), which are related through an element of the octahedral group that is 
not in the tetrahedral group. A sketch of that is given if Fig. 15. As can be 
seen there, adjacent vertices have independent tetrahedra, whereas all next- 
to-nearest neighbours correspond to the same tetrahedron and a relabelling of 
the poles. This relabelling does nothing for the Skyrme fields, and we are thus 
forced to identify these configurations as the same point. This tetrahedral in- 
variance means that the quantum mechanics is actually restricted to a quarter 
of the sphere, with slightly complex boundary conditions. The potential en- 
ergy has an even higher symmetry, however. Its fundamental domain, l/24th 
of the sphere, is mapped onto the whole sphere by the Oh group. Qualitatively 
one would expect that the lowest state has maximal symmetry. This is realized 
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if we make a symmetric combination of individual wave functions located in 
each minimum, leading to a cubic arrangement of wave functions. Note that 
this is not directly related to the pole positions; each minimum corresponds to 
a different tetrahedron. They can be arranged in two groups of four. Within 
each group the tetrahedra are related by a permutation of the poles (or an 
element of the tetrahedral group). The two groups are related by a 90° ro- 
tation, an element of the octahedral group. Since the pion field transforms 
non-trivial under all these transformations, we end up with a configuration 
explicitly symmetrised with respect to the full octahedral group. 



The tunnelling to the pretzel is somewhat harder to describe in the same 
terms, but that may not be necessary. We find that there are three - periodic 
- paths that describe tunnelling possibilities from one tetrahedron to one that 
is related by a reflection in a plane through the centre of the tetrahedron 
and parallel to two of the edges (an element of the octahedral and not of the 
tetrahedral group). Since these have a saddle point at lower energy than the 
paths through the doughnuts, one expects that they are even more important 
than those solutions. In essence these paths can be visualised as the edges of 
a cube, where each closed path is equivalent to one edge. 



Fig. 16. Labelling of the vertices of the cube. Each corresponds to a different real- 
ization of the tetrahedron 



If we label the different edges of the cube as in Fig. 16 we can estimate the 
spectrum using a simple Hiickel model, where we only include the nearest- 
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neighbour interaction diagonalising a simple 8-by-8 Hamiltonian matrix 



H 



( -E -X -X -X 

-X -E Q -X -X 

-X -E -X -X 

-X -X -E -X 

-X -E -X -X 

-X -X -E -X 

-X -X -£ -x 

yO -X -X -X -£ o y 



(45) 



where we assume all the eight configurations to be orthogonal. E is an esti- 
mate for the binding energy within a single well, and X is a mixing matrix 
element (whose sign will be positive). Since we have to identify even and odd 
tetrahedra, the only two eigenvectors of this matrix allowed are the symmetric 
and antisymmetric ones. The eigenvalues are — E ± 3X, and — E — 3X has 
eigenvector (1, 1, 1, 1, 1, 1, 1, 1). 



212 




100 



Fig. 17. The potential energy for a full period of motion along the circuit through 
the pretzels and the tetrahedron. 

Since this concerns a lower barrier, it is of greatest relevance to the large 
amplitude motion. The major issue here is that the harmonic oscillator length 
parameter for each of the wells is comparable to the length of the circuit for any 
reasonable value of Skyrme-model parameters. In this case the wave function 
in the circuit becomes approximately constant, and the energy is the average of 
the potential along the circuit. Since the contribution of the rotational kinetic 
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Fig. 18. The expectation value of the rotational energy for a full period of motion 
along the circuit through the pretzels and the tetrahedron. 

energy tends to flatten the potential even more, Fig. 18, this only reinforces 
this behaviour. If we just naively solve the one-dimensional problem on the 
circle, we find an energy of 209.5/ 7r /e + .00319/,,-e 3 , a considerable rise in 
energy. Furthermore, since the pretzel is a lot bigger than the tetrahedron, we 
find an enhanced value for the r.m.s. radius (by about 10%). 




(a) (b) 

Fig. 19. Two wave functions on the sphere, (a) x 2 y 2 + y 2 z 2 + z 2 x 2 and (b) 
xV + y 4 z A + z 4 x 4 + (x 2 y 2 + y 2 z 2 + z 2 x 2 )/2 

For the paths through the doughnut we cannot readily make such an estimate. 
Even though the harmonic oscillator length parameter is again close to the 
length of the path, which would imply a similar estimate as above, we cannot 
ignore the two-dimensional nature of the problem here. If we return to our 
VSEPR model, the simplest symmetric wave function, that has zeros at the 
infinities of the potential is x 2 y 2 + y 2 z 2 + z 2 x 2 (it is easier to write the wave 
function in terms of Cartesian tensors than as spherical harmonics). That 
wave function, as plotted in Fig. 19a, carries non-zero kinetic energy, and a 
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straightforward calculation shows that 

1 6 

(K) = -h\ (46) 

If we make the wave function more strongly centred along the tunnelling paths, 
as in Fig. 19b, the expectation value of K increases ( to 3.63fr 2 for this ex- 
ample). The best possible estimate for the contribution to the ground state 
energy from the two lowest modes is 

AE = -jp(K) + \J{V~ E )dl + i J E TOt dl. (47) 

Here R is the radius of the sphere associated with a tunnelling path of length 
L, L = 2 arctan(y^)i?. The integral is along one of the tunnelling paths. There 
are many objections one can make to this approximation (the real collective 
surface is not exactly a sphere, we ignore the zero-point motion orthogonal to 
the path, etc.) but we have been unable to come up with a reasonable alter- 
native. We estimate that (L 2 ) m 3.5h 2 , and calculate the remaining quantities 
from our results. 

Our final model is an independent quantisation of the five lowest non-zero 
modes, treating the remaining four in harmonic approximation. The potentials 
for the lowest five modes includes the difference of the local harmonic energy 
and the harmonic energy at the minimum, 

9 9 

IEMO)-iEMo). (48) 

i=5 i=5 

The rotational energy is treated in a similar way, i.e., we subtract the result 
at the minimum. The final result for the ground state energy is 

9 

E gs = tetrahedron + \ £ MO) + ^rot(O) + AE™ + 3AE&. (49) 

i=5 

Here we exhibit that we independent quantise modes three to five, and have 
an explicit two dimension problem for the lowest two modes. The various 
quantities in this expression are found to be (obtained in ways discussed above) 

E = 206.6—, 

e 

9 

i£M0) = 0.941Are, 

i=5 
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£ rot (0) = 0.00353^e 3 , 
A£ (12) = 3.8/ w /e - 0.037/^e + 0.0028/^e 3 , 
A£ (3) = 2.9/ 7r /e - 0.041/^e - 0.0003/^e 3 . 



(50) 



Taking all this together, and subtracting three times the B = 1 result, we find 



If we use our three parameter sets, we find values of 185.5 MeV for ANW, 
186.5 MeV for ANW and 211.6 MeV for R. These values are still high above 
the threshold, but please note that we have overestimated the effect of the 
lowest five modes by treating them as independent. Furthermore, some of 
the higher modes lead to a separated single hedgehog and a doughnut. Our 
experience from the B = 2 case suggests that in such cases the harmonic 
approximation may be relatively poor, and again leads to an overestimate. 
We have found a reduction of the energy, however, which is at least promising 
for a more mature approach to the problem. 



8 Conclusions and Outlook 

We have shown that our best treatment for the B = 2 problem leads to a state 
close to threshold (a few tens of MeV above or below) for a reasonable choice 
of parameters. At the moment this uses a enhanced version of the quantisation 
in the attractive manifold. What is really missing is the Hamiltonian in the 
full manifold, which is not readily accessible, even though one might consider 
cross-linking the three channels discussed in our paper, Ref. [1]. Using a dif- 
ferent numerical approach to study the modes leading away from the B = 2 
hedgehog, as done by Waindzoch and Wambach [20], may be another useful 
tool to get more information. It is unlikely that any of these methods will be 
able to completely cover the whole collective surface, however. 

A more pressing problem is the issue whether we can really construct a de- 
coupled 12-dimensional collective surface. Our results about level crossings in 
the attractive channel seem to suggest that it may not be possible to do so 
- or maybe we can, but we need to introduce a monopole geometric phase 
in our collective Hamiltonian, which should influence the dynamics and the 
rotational quantum numbers. This is an important issue in the quantisation 
of the B = 2 system. 

The B = 3 system apparently benefits a lot when the anharmonicities are 
taking seriously. In the present calculation we have not treated all of them. 
That is at least partly due to the fact that the JNR instantons are not the 



£ gs = -1.3^/6 + 0.781/^ 



- 0.007/ w e 3 . 



(51) 
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most general ones for this case. The fact that the lowest two modes separate 
in two bananas instead of three B = 1 Skyrmions as the monopole result from 
Ref. [17] seems to suggest should happen, is probably related to this deficiency. 

It is therefore necessary to see whether one can handle the AHDM instanton 
[21] in calculations of the type performed above. This is not trivial, since 
the relation between Skyrmion and AHDM data is rather indirect, and may 
just be too time consuming to implement. We are currently examining the 
feasibility of this approach. If we can we will be able to better address the 
anharmonicities for B = 3 and probably also see whether large amplitude 
fluctuations play a role for the B = 4 system. 
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A A detailed analysis of the modes of the tetrahedron. 

The instanton field is generated by four poles, arranged in an tetrahedron. All 
weights are equal. We choose 



where R is a parameter determining the size of the tetrahedron. 

Since all continuous symmetries are broken, we find 9 zero modes. In order 
of increasing frequency the non-zero modes can be interpreted as follows (we 
concentrate on how the pole positions change. Of course the weights change as 
well, but it is harder to classify modes that way.) Since we discuss fluctuations 
here, they have a natural identification as the gradients of a wave function. 
We shall identify those with the standard molecular orbitals found e.g., in 
Refs. [18,19]. We shall also apply the standard group representation labels 
corresponding to these wave functions. 



(Ai,Xx) 
(A 2 ,X 2 ) 
(A 3 ,X0 
(A4, X4) 



{l/4,-R,-R,-R,0), 



(l/4,-R,R,R,0), 
(1/4, R, -R,R,0), 
(1/4, R,R, -R,0), 



(A.l) 
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10-11 hoo = 0.318/^e, irrep=£. 

The two eigenvectors can be chosen to be 



S(Xi, Xi) — (0, -1, 1, 0, 0), 
5(X 2 ,X 2 ) = (0,-1, -1,0,0), 
5(A 3 ,X 3 ) = (0,1, -1,0,0), 

<5(A 4 ,X 4 ) = (0,1, 1,0,0), (A.2) 

and 

5(A 1 ,X 1 ) = (0, -1,-1,2,0), 
5(A 2 ,X 2 ) = (0,1, -1,-2,0), 
<5(A 3 ,X 3 ) = (0,-1, 1,-2,0), 

<5(A 4 ,X 4 ) = (0,1, 1,2,0). (A.3) 

They correspond to the gradient of the standard wave functions x 2 — y 2 and 
x 2 + y 2 — 2z 2 , respectively. The two modes, both invariant under a rotation 
around the x-axis, can be interpreted as follows. One mode corresponds 
to rotating the two lines (12) and (34) relative to one-another around the 
x-axis. The other correspond to stretching the lines (12) and (34), at the 
same time moving the lines closer. Both would be expected to go to a planar 
configuration of poles arranged in a square, the B = 3 doughnut. 
12-14 hto = 0.332/^e, irrep=T 2 . 

One representative of this mode, the one invariant under a 120° rotation 
about the 111 axis, has the form 

5(A 1 ,X 1 ) = (35A,1,1,1,0), 
<5(A 2 ,X 2 ) = (-5A,-1,0,0,0), 
S(X 3 ,X 3 ) = (-SX, 0,-1, 0,0), 

5(A 4 ,X 4 ) = (-<5A,0, 0,-1,0) (A.4) 

The wave function corresponding to the present mode is —\{yz + xz + xy), 
which has led to the T 2 assignment. 

In this mode we always keep the three points 2,3,4 in a shrinking equi- 
lateral triangle, at the same time moving pole 1 further and further away. 
This should probably connect to a well separated system consisting of the 
B = 2 doughnut and a single hedgehog. 
15-17 hu = 0.469/^e, irrep=T 2 . 

One representative of this mode, the one with eigenvalue 1 for a 120° rotation 
around the 111 axis, has the form 

<J(Ai,Xi) = (-3<JA,3,3,3,0), 
5(A 2 ,X 2 ) = (5A,-1,1,1,0), 
5(X 3 ,X 3 ) = (5X, 1,-1, 1,0), 

5(X 4 ,X 4 ) = (5X,1, 1,-1,0) (A.5) 
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The change in pole position is up to an overall translation in the (1, 1, 1) 
direction the same as found above. The important difference is that the size 
of 5X is much larger here (0.274), whereas the value found for the previous 
modes (.0383) is not incompatible with 0. 

Another interpretation is a rescaling of all poles, with a factor 1 + a for 
poles 2,3,4 and 1 — 3a for pole 1. 
18 ku> = 0.483/U, irrep=Ai. 

This is the breathing mode, as can easily be seen from the mode-vector. 
19-21 huj = 0.813/^e, irrep^Zl. 

One representative of this mode, the one with eigenvalue 1 for a 120° rotation 
around the 111 axis, has the form 

5(A 1 ,X 1 ) = (0,0,0,0,3), 
S{\ 2 , X 2 ) = (0, 0, .6278, -.6278, 1), 
<5(A 3 ,X 3 ) = (0, -.6278,0, .6278,1), 
5(A 4 , X A ) = (0, .6278, -.6278, 0, 1), 

8{c) = (1.178, 1.178, 1.178) (A.6) 

This is the only mode that mixes with iso-spin rotations. The change in 
the space components is consistent with the irrep T\ and a wave function 
x(y 2 — z 2 ) + y(z 2 — x 2 ) + z(x 2 — y 2 ). After a study of the change in baryon 
number associated with this mode we conclude that it probably corresponds 
to a twist, where we turn the baryon density in the 234 triangle, keeping 
the edges fixed at the point 1. At the same time a complicated grooming of 
the triangle relative to the pole 1 takes place (time-components). 
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